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1.  INTRODUCTION: 


The  objective  of  this  three  year  program  is  to  compute  the  center 
frequency  transmission  line  parameters  (namely  DELTERV  and  DELTERZ,  Ref.  1) 
for  electrodes  (including  buried  electrodes)  in  a  finite  length  single  or 
double  electrode  transducer  array.  These  parameters  are  to  be  calculated 
from  electrode  geometry  and  material  constants  taking  mass  and  elastic  loading, 
stored  energy,  piezoelectric  loading  and  external  circuits  into  effect. 

During  the  first  year  a  full  field  theory  was  developed  to  calculate 
all  of  the  effects  simultaneously.  During  the  present  year  a  different  approach 
was  adopted  in  order  to  obtain  transmission  line  model  parameters  that  are 
easily  expressed  in  terms  of  electrode  material  and  geometry.  The  scattering 
process  was  separated  into  three  categories  and  each  was  treated  individually 
using  perturbation  theory: 

*First  order  Mass  and  Elastic  loading 

*Stored  Energy  Effects  (Second  order  Mass  and  Elastic  loading) 

*Piezoelectric  loading. 

These  were  first  analysed  for  a  single  electrode  in  an  infinite  array. 

Section  2  describes  the  first-order  mass  and  elastic  loading  analysis  and 
Section  3  describes  the  analysis  of  acoustic  stored  energy  scattering.  The 
piezoelectric  scattering  is  analysed  in  Section  4.  This  analysis,  however, 
does  not  allow  for  interconnections  between  electrodes  and  so  is  not  used 
directly  for  interdigital  transducers.  In  Section  5  an  analytical  expression 
for  the  frequency  response  of  a  single  electrode  in  an  infinite  grounded  array 
is  determined.  The  results  of  Sections  4  and  5  have  been  combined  to  get  a 
frequency  dependent  scatter  matrix  for  a  single  electrode  in  an  infinite  array 
excited  with  arbitrary  voltages  on  the  electrodes.  This  will  be  described  in 
the  final  report.  The  infinite  array  analysis  for  the  contract  requires  the 
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scattering  parameters  to  be  evaluated  only  at  harmonics;  however,  the 
finite  array  analysis  presently  under  development  requires  the  scattering 
parameters  to  be  evaluated  at  all  frequencies. 


r 


All  of  these  scattering  contributions  are  combined  to  provide  a 
composite  mixed  unit  scatter  matrix  description  for  the  electrode.  The  des¬ 
cription  is  left  in  the  scatter  matrix  form  because  the  transmission  line 
model  for  odd  and  even  harmonic  scattering  is  different.  This  scatter  matrix 
is  easily  converted  to  the  DELTERV  and  DELTERZ  for  the  odd  harmonic  trans¬ 
mission  line  parameters  using  the  equations  derived  in  Appendix  I. 


i 


2.  FIRST  ORDER  ELASTIC  LOADING: 


A.  Introducti  .1: 

The  problem  of  surface  wave  reflection  at  thin  strip  overlays 

occurs  frequently  in  signal  processing  devices.  Grooved  array  reflectors 

employ  strip  overlays  of  the  same  material  as  the  substrate,  while  metallic 

transducers  incorporate  strips  of  a  material  different  from  the  substrate. 

(2) 

In  his  book,  Auld  has  described  a  method  for  determining  the  first-order 

reflection  coefficient  of  a  very  thin  strip  overlay.  In  this  method  the 

normal  and  tangential  stresses,  TS  ,  TS  ,  TS  (Fig.  lb),  generated  at  the 

xy  yy  zy 

interface  by  the  overlay  are  used  as  source  terms  in  the  normal -mode 
equation  for  the  generation  of  the  reflected  wave.  The  results,  however, 
when  applied  to  practical  cases  such  as  grooves  in  ST-Z  quartz  or  aluminum 
on  ST-X  quartz,  yield  a  value  of  the  reflection  coefficient  that  is  too 
high  compared  to  experimental  values. 

s  s 

In  Ref.  2  ,  the  interface  stresses  under  the  overlay  T  ,  T  , 

xy  yy 

s 

T  are  assumed  to  be  the  same  as  those  due  to  a  uniform  overlay  extending 
all  along  the  propagation  path.  In  this  section  it  is  shown  from  the  first 
order  equations  that,  in  addition  to  this  uniform  component,  a  pair  of  delta 
function  shear  stresses  is  generated  at  the  edges  of  the  strip  due  to  its 

finite  extent.  The  existence  of  these  stresses  has  also  been  shown  recently 

(8 ) 

by  a  rigorous  perturbation  scheme.  By  incorporating  the  additional  source 

terms  into  the  normal  mode  equations,  the  calculated  reflection  coefficients 
for  grooves  in  ST-X  quartz  and  for  aluminum  on  ST-X  quartz  are  found  to  be 
in  close  agreement  with  observed  experimental  values.  The  reflection  coeffi¬ 
cient  calculated  for  grooves  in  Y-Z  lithium  niobate  also  agrees  with  the 
experimental  value;  however,  in  this  paper  any  piezoelectric  effect  due  to 
the  overlay  is  neglected.  This  effect  is  negligible  in  quartz  but  not 


necessarily  so  in  lithium  niobate. 


The  present  formulation  thus  affords  a  simple  technique  for  com¬ 
puting  the  first-order  reflection  coefficient  directly  from  material  para¬ 


meters  . 


B.  Formulation: 

The  incident  surface  wave  (Fig.  la)  is  assumed  to  travel  in  the 

positive  z-direction  in  a  half-space  with  a  free  surface  so  that  the  surface 

stresses  T  ,  T  ,  T  are  zero.  The  wave  is  uniform  in  the  x-direction 
xy  yy  zy 

and  travels  with  a  velocity  c^  such  that  the  wave-number  {5  is  given  by 
P  =  (u/c 

s 

where  u>  is  the  radian  frequency.  A  thin  strip  of  height  h  extends  along 

A  2L 

the  propagation  path  from  z  =  -  —  to  z  =  +  -  (Fig.  la).  Due  to  the 

s  s  s 

presence  of  the  strip,  there  are  interface  stresses  T  ,  T  ,  T  (Fig.  lb) 

xy  yy  zy 

(2) 

which  are  the  source  terms  in  the  normal  mode  equations  describing  the 

change  in  the  incident  forward  wave  and  the  growth  of  a  reverse  wave. 

Let  A+  and  A  be  the  amplitudes  of  the  forward  and  reverse  waves. 

The  surface  wave  in  general  has  three  particle  velocity  components  v  ,  v  , 

x  y 

v  whose  values  at  the  surface  for  the  forward  wave  are  written  as, 
z 


(v+  e'jez) 

X 


(v+  e'jPz) 

y 


(la) 

(lb) 

(lc) 


R  R  R 

The  particle  velocities  for  the  reverse  wave  v  ,  v  ,  v  are  the 

x  y  z 

(3) 

negative  complex  conjugates  of  those  of  the  forward  wave. 


(2b) 


R  +  -jBz  * 

v  =  -A  (v  e  ) 

y  y 

R  +  -jPz  * 

v  =  -A  (v  e  ) 

z  z 


(2c) 


The  normal  mode  equations  are: 


1 

dA+  = 

1 

'  * 

A+ 

dz 

|A+|2.4Pa 

V 

L  x 

_1_ 

dA~  = 

1 

~  R* 

A" 

dz 

1 

V 

L  x 

S  *  s  *  „s 

T  +  v  T  +  v  T 
xy  y  yy  z  zy_ 


TS  +  vR*  TS  +  vR*  TS 
xy  vy  yy  z  zyJ 


(3) 


(4) 


where  P  is  the  total  power  per  unit  beamwidth  in  a  unit  amplitude  surface 

3L 

wave. 

Equation  (3)  describes  the  change  in  velocity  of  the  forward  wave. 
The  first  order  theory  to  be  described  here  yields  the  same  velocity  change 
as  in  Ref.  2;  so  this  is  not  discussed  further  in  this  paper.  Equation  (4) 
describes  the  growth  of  the  reflected  wave.  Using  Equations  (1)  and  (2), 
it  is  written  as, 

s 


1  dA~ 

.+  dz 

A 

[■i  (¥) 

X 

where 

/  +  J 

(V 

j2Pz 


4P 


(3) 


(6a) 


,  -N2 
4P 


(6b) 


,  -s2 

(v  ) 
z 

4P 


(6c) 


3i  fX  t 

Integrating  Equation  (4)  from  z  =  -  —  to  +  ^  and  assuming  A  to  be  constant 


over  this  region  we  obtain  the  reflection  coefficient,  r. 


6 


where 


xx  ■  f  {¥)  *-j2te  - 


(8a) 


_a  x 
2 


‘y  ■  ?  (>)  *'32te  <■ 
3  a  y 

'2 


(8b) 


i.  ‘  f  &)  <■ 

3  z 
'2 


(8c) 


C.  Interface  Stresses  Generated  by  Strip  Loading: 

To  obtain  1,1,  and  I  ,  we  have  to  determine  the  induced  stresses 
x  y*  z* 

s  s  s 

T  ,  T  ,  and  T  at  the  interface  between  the  strip  and  the  substrate.  In 
xy  yy  zy 

the  first-order  approximation,  the  particle  velocities  throughout  the  strip 

can  be  assumed  to  be  the  same  as  those  at  the  surface  of  the  substrate  due 

,  (2) 

to  the  incident  wave.  The  interface  stresses  are  then  readily  determined 
by  considering  the  forces  acting  on  a  differential  element  dx*  z  of  the 
strip  of  height  h.  From  the  x-directed  forces  (Fig.  2a), 


ar' 

T9  =  juip'hv  -  h. 
xy  J  r  x  dz 


(9a) 


From  the  y-directed  forces  (Fig.  2b), 


T  =  jrnp'hv 

yy  y 


(9b) 


From  the  z-directed  forces  (Fig.  2c), 


s  ST ' 

T  =  jrnp'hv  -  h.  zz 
Zy  2  ~ 

where  p'  is  the  mass-density  of  the  strip. 
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(9c) 


Fig.  2 


Traction  rree 
“■i  ,  Too  Surface 

Jx  <■' 


Ty'2  •  nax 


(b) 


First-order  forces  acting  on  a  differential  element  of  the  strip  in 

(a)  x-direction 

(b)  y-direction 

(c)  z-direction. 
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T'  has  been  neglected  In  Equation  (9b)  because  its  first  order 
yz 

value  Is  zero  In  the  strip.  The  primes  are  used  on  T'  T’  and  T'  to 

xz’  yz  zz 

denote  that  these  are  fields  In  the  strip  rather  than  In  the  substrate. 

These  results  have  also  been  derived  in  Ref.  2  using  a  rigorous  Taylor's 

series  expansion  for  the  fields  In  the  strip. 

In  Equations  (9a)  and  (9c)  T'  and  T'  are  the  stresses  in  the 

xz  zz 

strip  which  is  assumed  to  have  the  same  particle  velocities  v  and  v  as 

X  z 

the  substrate.  We  may  write, 


T 1 
xz 


=  -  &  a'  v 

Ctf  X  X 


(10a) 


i 


T '  =  -  j-  v 

ZZ  U)  Z  Z 


(10b) 


where  or'  and  or'  are  constants  with  the  dimensions  of  stiffness  (nt/m  ). 
x  z 

These  are  effective  stiffness  coefficients  that  depend  on  the  material  of 
the  strip.  In  a  strip  made  of  an  isotropic  material  with  Lame'  constants 
X'  and  p' ,  a£,  are  readily  obtained  from  the  plane  strain  equations 
assuming  the  first  order  T'  to  be  zero. 

yy 


or'  =  p.' 

X 

( 1  la) 

(lib) 

X'  +  2  pi* 


For  an  anisotropic  strip  equations  similar  to  (11)  are  derived  but  are  more 
(21 

complicated.  '  '  However,  for  an  anisotropic  strip  of  the  same  material  as 

the  substrate,  T'  and  T'  are  the  same  as  the  surface  values  of  these 
zz  zx 

quantities  in  the  incident  unperturbed  surface  wave. 


T'  =  T 
zz  zz 


(12a) 


T'  =  T  (12b) 

xz  xz 


where  T  and  T  are  the  surface  stress  fields  ot  the  incident  surface  wave. 


It  will  be  noticed  that  the  derivatives  of  T'  and  T'  with  respect 

xz  zz 

to  z  appear  in  Equations  (9a)  and  (9c).  These  stresses  have  a  uniform  value 
(multiplied  by  the  factor  exp[-jf5z])  over  the  entire  length  of  the  strip. 

Si  cl 

However,  at  the  end  faces  (z  =  -  ^  and  z  =  +~)  these  stresses  must  be  zero. 
This  leads  to  the  excitation  of  evanescent  modes  around  each  edge  which  give 

2  (g) 

rise  to  the  second  order  (~  Oh)  )  stored  energy  effects.  In  the  present 

first-order  approximation  we  will  neglect  these  evanescent  modes  and  assume 

that  the  stresses  T'  and  T'  remain  uniform  over  the  length  of  the  strip 
xz  zz 

and  abruptly  go  to  zero  at  the  edges.  Because  of  this  abrupt  change,  the 
derivative  terms  in  Equations  (9a)  and  (9c)  give  rise  to  delta  function 
stresses  at  the  two  edges.  We  thus  have  from  Equations  (9)  and  (10), 


T 1  =  KP'Cs  -  “x)  +  jf  (5(Z  +  2)  -  6(2  -  fOj 


(13a) 


yy  _  i(gh) 


[■•■a 


(13b) 


[(*•■;  -«i)  (•<■♦!> -•<■-!>)] 


(13c) 


where  c  =  uu/p  is  the  velocity  of  the  surface  wave, 
s 

Fig.  3  shows  the  spatial  distribution  of  the  interface  stresses. 

s  s 

The  delta  function  components  of  T^  and  T  at  the  edges  of  the  strip  were 
not  included  in  Ref.  1.  Inserting  the  computed  stresses  from  Equations  (13) 
into  Equations  (8)  we  get. 


imi  (p.c2  +  a.)  -iafia 

C  K  s  X  P 


v  -  -4P  ^ 


(14a) 


(14b) 


(14c) 


I 

z 


=  Jim 

c 

s 


<p'% -  «;> 


Comparing  Equations  (13)  and  (14),  it  will  be  noticed  that  the  effect  of 
adding  the  delta  function  terms  is  to  change  the  sign  of  a £  and  a '  in 
Equations  (13a)  and  (13c).  Physically,  this  means  a  change  in  the  sign  of 


the  component  of  reflection  die  to  stiffness  since  a^,  represent  stiffness 
while  p’  represents  mass. 


(b) 


Fig.  3.  Spatial  distribution  of  (a)  x-directed,  (b)  y-directed,  and 


(c)  z-directed  interface  stresses. 


Combining  Equations  (14)  with  Equation  (7),  we  have, 


Fn'c2 

(Zl  +  Z2  +  Z3\ 

4. 

/Of'z  +  a'zoX 

/  x  1  z 

Lp  cs 

\  0)  / 

T 

VO)/ 

D,  Application  to  Specific  Cases: 

In  this  section  Equation  (15)  is  applied  to  obtain  the  first-order 

reflection  coefficients  per  strip  for:  (1)  aluminum  on  ST-X  quartz  and  (2) 

Grooves  in  ST-X  quartz.  In  each  case  the  computed  value  closely  matches 

the  experimental  value  reported  in  Ref.  4. 

The  values  of  relevant  field  quantities  at  the  surface  for  the 

unperturbed  surface  wave  in  ST-X  quartz  are  given  below  for  a  total  power 

12 

per  unit  beamwidth  P^  =  10  ui  Watts/m.  These  values  were  computed  from  the 
field  equations  using  the  method  in  Ref.  5. 


v+  =  .40  uj  <  -90°  m./Sec. 

x 

v+  =  4.24  tu  <  +90°  m./Sec. 

y 

v+  =  2.79  uj  <  0  m./Sec. 

z 

T  =  2.42xlOU  g  <  180°  nt./nn 

zz 

10  ? 

T  =  2.66x10  P<  90°  nt./rn. 

zx 


The  mass-density  and  surface  wave  velocity  for  quartz  are, 

p  •=  2651  kg/m? 

c  =  3 158 . 5  m. /Sec . 
s 

We  thus  have, 

Z\  /  2. 

—  =  -4x10  m./nt. 
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1 


Z2  -12  2 
■—  =  -4.5x10  mf/nt. 

Z1  -1?  2 
—  =  +1.95x10  m./nt. 

0} 

Using  these  values  in  Equation  (15)  we  get, 

.  r  i  of'  -i 

r  =  j(sinpa)  .  ^  I -.43  +  .0066  +  .32  J  (16) 

Pcs  PCs 

This  equation  gives  the  reflection  coefficient  of  a  strip  of  any  material 
on  ST-X  quartz.  The  strip  material  is  characterized  by  p',  a^,  and  a^,  p 
being  the  mass -density  of  quartz. 

(1)  Aluminum  on  ST-X  quartz: 

Aluminum  has  Lame'  constants  X1,  p.'  and  mass-density  p'  equal  to, 

X'  =  6.1xl010  nt./nu 

10  2 
p,'  =  2.5x10  nt./m. 

p'  =  2695  kg./m? 

Using  Equations  (11a)  and  (lib), 

10  2 

at'  =  7.75x10  nt./m. 
z 

10  2 

a'  =  2.50x10  nt./m. 
x 

We  thus  have, 

p'/p  =  1.017 

a'/pc2  =  0.945 
x  's 

or'/pc2  =  2.930 

Z  r  3 
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Hence  from  Equation  (16), 


r  =  j(singa)  .  (.507  7) 

For  a  reflector  array  at  its  fundamental  stopband  with  equal  electrode  and 
gap  regions,  singa  =  1,  so  that  the  reflectivity  per  strip  is  .507 

(2)  Grooves  in  ST-X  quartz: 

In  this  case  the  overlay  is  of  the  same  material  as  the  substrate 

as  discussed  in  Section  C,  T'  and  T'  in  the  strip  may  be  assumed  to  be 

zz  zx 

the  same  as  their  :urface  values  for  the  unperturbed  waves  which  are  listed 
at  the  beginning  of  this  section.  Using  these  values, 

T 

or  =  -  -r  -  =  8.67x10  nt./m. 

z  g  v 

z 

,  Txz  ,  ,,  ln10  ,  2 

of  =  -  —  -  =  6.65x10  nt./m. 

X  g  v 

X 

We  thus  have, 

p'/p  =  1 

of’/pc2  =  2.514 

X  s 

Qf'/pc2  =  3.278 
z  s 

Hence  from  Equation  (16), 

r  =  j (singa)  (.636  ^) 

A 

A  similar  calculation  with  Y-Z  lithium  niobate  as  the  substrate  material 
yields, 

,  p  ,  or  •  1 

r  =  j  (singa)  .  -  .317  +  .272  — ^  J  (17) 

P 


14 


where  p'  and  a '  describe  the  strip  material  and  p  is  the  mass -density  of 
2  10  2 

lithium  niobate  (Pcg  =  5.72x10  nt./m.).  does  not  appear  in  this  equation 

because  v  =0  for  a  surface  in  Y-Z  lithium  niobate.  For  grooves  in  Y-Z 
x 

lithium  niobate  we  have, 

a'  =  -  f  —  =  2.092xl0U  nt./m? 
z  p  v 

z 

so  that, 

p'/p  =  1 

of*  /  pc2  =  3.657 
z  v  s 

Equation  (14)  then  yields, 

r  =  j(singa)  .  (.678 

A. 

This  result  too  agrees  with  experimental  data.  However,  as  noted  earlier, 
we  have  neglected  the  electrical  term  in  the  normal  mode  equation. 

E.  Conclusions: 

In  this  section  the  first-order  reflection  coefficient  of  strip 
overlays  on  anisotropic  substrates  is  obtained  from  the  well-known  normal  mode 
theory.  An  analytical  expression  for  the  reflection  coefficient  is  derived 
for  arbitrary  anisotropic  strip  and  substrate  materials.  The  computed  values 
agree  well  with  experimental  values  for  grooves  in  ST-X  quartz  and  Y-Z 
lithium  niobate  and  for  aluminum  on  ST-X  quartz. 
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3.  STORK!)  ENERGY  SCATTERING  : 


A.  Introduction: 

The  propagation  of  surface  acoustic  waves  through  a  periodic  array 
of  grooves  or  thin  strip  overlays  has  been  modeled^'*  as  a  repetitively  mis¬ 
matched  transmission  line  with  reactive  energy  storage  elements  ( j B )  at  each 
discontinuity  (Fig.  4).  The  mismatch  parameter  AZ/Z  has  been  related  to  the 
elastic  constants  of  the  substrate  and  strip  material  using  various  theore¬ 
tical  approaches . ^ ^ ^  However,  no  theoretical  analysis  relating  the 
energy  storage  elements  to  material  constants  has  yet  been  described. 

The  energy  storage  element  (jB)  is  physically  attributed  to  the 

2 

excitation  of  evanescent  modes.  It  is  a  quadratic  effect  of  order  ~  (h/\)  , 
where  h  is  the  height  of  the  strips;  this  is  in  contrast  to  AZ/Z  which  is  a 
linear  effect  of  order  ~  (h/\).  This  energy  storage  produces  two  specific 
effects  on  the  dispersion  curve  which  have  been  verified  experimentally:  (1) 
it  leads  to  a  shift  in  the  center  of  the  fundamental  stop-band,  and  (2)  it 
predicts  a  stop-band  at  the  second  harmonic. 

In  this  section  we  will  describe  a  model  using  a  first-order  pertur¬ 
bation  analysis  of  the  strip  overlay  that  predicts  these  effects  and  shows 
their  physical  origin.  Piezoelectric  effects  are  neglected  in  this  paper; 
only  mechanical  effects  are  considered.  However,  anisotropy  is  accounted 
for.  The  results  for  a  few  practical  cases  involving  aluminum  electrodes 
on  ST-X  quartz  and  grooves  in  ST-X  quartz  are  in  fair  agreement  with  re¬ 
ported  experimental  values. 

B.  Secular  Equation: 

The  field  quantities  for  a  wave  propagating  in  a  periodic  structure 
satisfy  Floquet's  Theorem. ^  ^  The  particle  displacements  and  the  stresses  at 


1  < 

1  i  1 

1  Surface  Wave  |  j 

Tj-  a/p  J 

! 

[ 

™  J l 

L  -'p/zJ 
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m 


K-O*^ 

—  P- 


•n  a  periodic  array 
,  Sur£ace  wave  propagation  m 
Fig.  bur 

(a)  Physical  Structure 

(b)  Equivalent  Circuit 
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the  interface  between  the  layer  of  strips  and  the  substrate  may  by 


written  as. 


+  CD 


IK  =  ^  exp(-jP^z) 

n=  -co 
f  oo 

\  (n)  .  ,a(n)  . 

Ti  =  2^  i\  exp(-jp  z) 


n= 


where,  -  0  +  n.2-rc/p 


(18a) 


(18b) 


(19) 


p  being  the  period  of  the  strip  array.  Here  represents  the  stress  at  the 
interface.  The  second  subscript  is  implied  to  be  y  and  is  dropped  for  con¬ 
venience.  p  is  the  surface  wave  wavenumber  and  represents  the  n'"*1 

space  harmonic. 

The  relationship  between  r  and  U  required  by  the  periodic  array  of 
strips  is  determined  using  a  first-order  perturbation  analysis;  this  analysis 
is  basically  similar  to  that  described  by  Auld  in  his  book^1^  with  a  modifi¬ 
cation  to  handle  edge  effects . ^  Using  this  analysis,  a  matrix  [L],  charac¬ 
teristic  of  the  strip  material  and  geometry  (Fig.  5),  is  evaluated  in 
Appendix  II  such  that. 


4- co 


T(m)  =  CPh)  .  )  )  L(m’n)  U(n) 

1  n=-~  j=x  ij  j 


l  =  x,y,z 
m  =  -  a>,  oo 


(20a) 


where  h  is  the  height  of  the  strips.  A  first-order  perturbation  analysis 
yields  stresses  that  depend  linearly  on  (Ph).  The  strip  array  being  periodic 
couples  together  different  space  harmonics.  A  particle  displacement  at  one 
spatial  harmonic,  n  excites  stresses  at  another  spatial  harmonic,  m. 

This  is  in  contrast  to  the  uniform  substrate  which  does  not  couple 
different  space  harmonics.  A  matrix  { S ] ,  characteristic  of  the  substrate 
material  (Fig.  b)  ,  is  evaluated  in  Appendix  III  such  that, 
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Fig.  5.  Stress-displacement  relation  for  a  semi-infinite  substrate  and 
for  a  periodic  array  of  strips. 
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r  - 1  i  *Jr!  ujn>- 

ra  —  _  03  1  =V  J 


n=  -®  j=x 


l  =  x,y,z 


m  = 


(20b) 


Since  the  uniform  substrate  does  not  couple  different  space  harmonics,  we 


s(m,n)  =  0  for  m  ^  n  (21) 

J-J 

The  [Lj  and  [Sj  matrices  thus  describe  the  elastic  properties  of  the  periodic 
layer  of  strips  and  the  uniform  substrate  respectively.  The  [S]  matrix  repre¬ 
sents  an  exact  analysis  of  the  semi-infinite  substrate,  while  the  [L]  matrix 
is  based  on  a  first-order  perturbation  analysis  of  strips. 

The  stresses  at  the  interface  are  equal  and  opposite  for  the  sub¬ 


strate  and  the  strips,  so  that  equations  (20a)  and  (20b)  are  combined  to  yield. 


-4-oo  7 

\  \  Pc (m, n) 


n=  J=x 


^  [g(m,n)  +  (ph)>L(m,n)j  y(n)  =  Q 

1  — v  — *  f  J  J  1 


i  =  x,y,z 


m  =  -<*>,<*> 


(22a) 


In  matrix  form, 


[[S]  +  (ph).[L]J  {U}  =  {0} 


(22b) 


where  [u]  is  a  3M  x  1  column  vector  and  [S]  and  [L]  are  3M  x  3M  matrices,  if 
M  is  the  number  of  space  harmonics  being  considered.  The  factor  3  arises 
because  of  the  three  dimensions,  x,  y,  and  z. 

Equation  (22b)  is  the  secular  equation  for  the  infinite  set  of  space 
harmonic  amplitudes.  The  elements  of  [Sj  and  [ L]  are  functions  of  P  and  u>  so 
that  in  principle,  one  can  solve  for  the  dispersion  relation  P(u>)  in  any  degree 
of  truncation  approximation  by  setting. 


det  [SLJ  =  0 


(23a) 


where 


tSL]  =  [S]  +  (ph)  .  [L] 


(23b) 
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In  the  limit  h  -•  0,  the  strips  vanish  and  there  is  no  coupling  between  the 
space  harmonics.  Equation  (23a)  then  becomes  det  [S]  =  0  and  this  yields 
the  dispersion  curve  for  surface  waves  in  a  substrate  with  a  free  surface. 


C.  Classification  of  Space  Harmonics: 

Before  we  discuss  the  solutions  to  equation  (23a),  it  is  useful 
to  give  a  physical  interpretation  of  the  different  space  harmonics.  We 
note  that. 


P 


(24a) 


is  the  wavenumber  of  the  forward  propagating  surface  wave.  If  P  is  at  one 
of  the  stopbands,  we  may  write, 

P  =  N  .  n/p  (24b) 

where  N  is  an  integer.  N  =  1  represents  the  fundamental  stopband,  n  =  2 
represents  the  second  harmonic  stopband  and  so  on.  We  then  have  from 
equation  (19), 

P('N)  =  -P  (24c) 


so  that  the  space  harmonic  (-N)  represents  the  backward  propagating  surface 
wave.  As  may  be  expected,  at  the  stopbands  the  forward  and  backward  waves 
are  coupled  together. 

Now,  the  space  harmonics  with  -N  <  n  <  0  have  |p*‘n^|  <  |p|,  and 

thus  couple  to  propagating  bulk  modes  that  carry  power  away  from  the  surface. 

The  space  harmonics  with  n  <  -N  or  n  >  0,  on  the  other  hand  have  } |  >  |p|. 

They  couple  to  evanescent  modes  that  store  power  reactively  near  the  surface. 

To  summarize: 

(1)  n  =  0  represents  the  forward  wave. 

(2)  n  =  -N  represents  the  backward  wave. 
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(3)  -N  <  n  <  0  represents  propagating  bulk  modes  that  cause  energy  loss. 

(4)  n  <  -N  and  n  >  0  represent  evanescent  modes  that  cause  energy  storage. 

This  is  illustrated  in  Fig.  6  for  N  =  3.  It  will  be  noted  that  at  the  funda¬ 
mental  stopband  N  =  1  so  that  there  are  no  space-harmonics  in  category  3.  It 

is  for  this  reason,  we  believe,  that  bulk  power  losses  are  much  larger  for 

resonators  operated  at  the  higher  harmoincs  than  for  those  operated  at  the 

t  undamenta 1  stopband . 
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Fig.  6.  Spatial  harmonics  excited  by  an  incident  wave  in  a  periodic  array 
at  the  third  harmonic  stopband  (N  =  3). 


D.  Coupled-Mode  Approximation: 

In  solving  Equation  (23a)  for  the  secular  determinant  it  is  common 
to  use  the  coupled-mode  approximation  near  a  stopband.  At  the  stopband 
the  space  harmonic  amplitudes  for  n  =  0  and  n  =  -N  (representing  the  forward 
and  backward  space  harmonics)  are  the  largest  and  only  these  two  are  retained 
We  then  write, 


det 


SL 


(0,-N) 


Ls 


l(-h,°)  sl(-n,-n* 


=  0 


(25) 


Here  only  two  space  harmonics  are  retained  and  the  matrix  is  really  6x6  since 
3  dimensions  are  associated  with  each  space  harmonic.  Thus  Sl/^’  ^ 
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Si/  N,<^')  and  Si/  N>  represent  (3x3)  matrices,  but  for  clarity  these  are 
written  as  single  terms. 

It  will  be  noted  that  the  off-diagonal  terms  in  Equation  (25)  come 
from  the  periodic  strip  array.  The  substrate  matrix  [S]  is  zero  between 
different  space  harmonics.  Thus  the  off-diagonal  terms  are  all  of  order 
h/X.  For  example, 


sl(0’-n) 


Oh) 


l(°,-N) 


(26) 


since 

s<0--N>  -  0 

The  dispersion  curve  can  be  obtained  from  Equation  (25).  Equivalently  the 

(12) 

normal-mode  theory  can  be  used.  Either  approach  yields  the  mismatch  re¬ 
flection  represented  by  AZ/Z,  and  the  results  agree  with  other  computa- 
tions^’^’^  and  with  experiment.  But  this  coupled  mode  approximation  does 

not  predict  the  "stored  energy"  effect;  there  is  no  shift  in  the  centers 

2 

of  the  stopbands  of  order  (h/X)  .  At  the  second -harmonic  (N  =  2),  for  equal 
strips  and  gaps  it  is  found  that. 


^(0,-2) 


l(-2,0) 


0 


(27) 


which  means  that  no  stopband  is  predicted  at  the  second  harmonic  since  for¬ 
ward  and  backward  waves  are  effectively  decoupled  in  this  approximation. 

It  will  be  noted  that  carrying  the  perturbation  to  second-order  in 
the  analysis  of  the  strip-array  (that  is,  in  calculating  the  matrix  [ L] )  does 
not  yield  the  stored  energy  effect;  there  is  still  no  coupling  between  the 
forward  and  backward  waves  at  the  second  harmonic. 

To  summarize,  the  coupled-mode  approximation  gives  us  the  AZ/Z  in 
the  transmission  line  model,  but  not  jB  (Fig.  4).  Carrying  the  perturbation 
analysis  of  the  strips  to  higher  orders  only  improves  AZ/Z  to  higher  orders. 
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E.  Beyond  the  Coupled-Mode  Approximation: 

To  obtain  the  "stored  energy"  effect  one  needs  to  go  beyond  the 
coupled  mode  approximation.  To  illustrate  the  point,  let  us  include  one 
space-harmonic,  n  in  addition  to  0  and  -N  representing  the  forward  and 
backward  harmonics.  Equation  (22b)  is  written  as. 


'SL<°’0) 

sl(o,-m> 

o 

D 

- — •• 

_ 1 

u(0> 

sl(-n’0) 

sl(-n,-n) 

U(-N> 

=  0 

(28) 

SI (n,0) 

SL(n,-N) 

SL(n,n) 

u(n) 

From  these  equations  is  eliminated  to  yield, 

'SL.(0,0)  sl,(0,-N)1  jj  ( o ) 

=  0  (29) 

sl,(-n,°)  sl,(-n,-n)  u(-h) 


where 

Sl/0'0' 

=  SL<°'0) 

-  SL(°*n) 

Si/0'’"' 

=  sl(0’-n) 

-  SL(°’n) 

SL'<-,:'0) 

=  sl(-n’0) 

-  SL(-N'n) 

=  sl(_n’"n) 

-  SL(-N’n) 

[SL(n’n)]-1 

SL(n’0) 

(30a) 

[  SL(n>  n) ] -1 

SL(n’-N) 

(30b) 

[SL(n,n)  ]  _1 

SL(n’0) 

(30c) 

( SL(n,n)  ] _1 

SL(n’-N) 

(30d) 

Equation  (29)  is  of  the  same  form  as  obtained  from  the  coupled-mode  approxi¬ 
mation  but  the  matrix  elements  are  modified  as  give  by  Equations  (30).  The 
dispersion  curve  is  obtained  by  setting, 

'Sh,(0>0)  SI.,(°’'N)" 

det  =  0  (31) 

sl,(-n,°)  sl.(-n,-n) 

but  now  the  effect  of  the  space  harmonic  n  is  included  by  using  the 
corrected  matrix  elements  SI/  instead  of  SL. 
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It  will  be  noted  that  the  correction  terms  in  liquations  (30)  are 
2 

each  of  order  (h/X)  since  they  contain  the  product  of  two  off-diagonal  terms, 
each  of  which  is  of  order  h/X  as  discussed  previously.  Moreover,  at  the 
second -harmonic  stopband  SL'^>  ^  o  (though  SL^’  ^  =  0)  because  of  the 

extra  term  in  Equation  (30b)  arising  from  the  additional  space  harmonic 
considered.  The  dispersion  curve  now  shows  a  shift  in  the  center  of  the 

stopband  (through  SL'^°’^)  and  also  has  a  stopband  at  the  second-harmonic 

(0-2)  2 
(through  SL'  ’  ).  These  effects  are  of  order  (h/X)  and  arise  from  the 

inclusion  of  the  extra  space  harmonic  representing  a  propagating  or  an 

evanescent  mode  as  discussed  previously.  The  corrections  in  Equations  (30) 

are  readily  generalized  to  include  as  many  other  space-harmonics  as  desired: 


SL' 

(0,0) 

SL(°’’ 

0) 

-  z 

SL(°’n) 

G(n,m) 

SL(m’0) 

n,ra 

SL' 

(0,-N)  = 

SL<°’ 

-N) 

-  Z 

SL(°»n) 

G(n,m) 

SL(m,-N) 

n,m 

and  similarly 

for  SL' 

(-N,0) 

and 

SL '  ^  ” 

where, 

__  -1 

[G]  =  [SL]  (33) 

[SL]  being  the  matrix  [SL]  with  the  rows  and  columns  corresponding  to  0  and 
-N  deleted.  The  summation  over  n,m  can  be  carried  out  over  as  many  space- 
harmonics  as  desired  and  the  total  "stored  energy"  effect  is  thus  obtained. 

Equation  (31)  can  then  be  solved,  in  principle,  to  obtain  the 

(12) 

dispersion  curve.  Instead  we  will  use  the  normal  mode  theory  to  calculate 
the  phase-shift  and  the  reflection  per  periodic  section.  This  is  simpler  to 
implement  than  the  direct  solution  of  the  secular  equation. 

F.  Normal  Mode  Approach: 

In  the  secular  determinant  approach,  we  determine  the  appropriate 


u)  for  a  given  g  so  that  equation  (31)  is  satisfied.  A  non-trival  solution 


U(0) ,  U('N)  is  then  obtained  that  satisfies  Equation  (29).  Using  the  normal 


mode  theory  the  problem  is  approached  differently.  To  start  with  we  assume 


that  we  have  an  unpe rturbed  surface  wave  of  amplitude  A  (Fig.  7).  This 
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Fig.  7.  Surface  wave  scattering  from  a  single  strip  in  a  periodic  array. 


wave  is  a  normal  mode  of  the  semi-infinite  medium  so  that  u>  and  g  are  re¬ 
lated  by, 

id  =  gcs  (34) 

where  cs  is  the  velocity  of  the  surface  wave  on  a  free  surface.  It  will  be 
noted  that  the  dispersion  relation  (34)  is  a  solution  of 

det  IS]  =  0  (35) 


rather  than  a  solution  of  equation  (31). 

Since  this  wave  is  not  a  normal  mode  of  the  periodic  structure, 
Equation  (29)  is  not  satisfied  and  a  residual  stress  distribution  is  set  up. 
If  the  incident  forward  propagating  normal  mode  has  a  particle  displacement 

u  (u  being  a  vector  with  three  components)  then  the  space  harmonics 
(-N) 

and  t  of  the  residual  stress  are  given  by. 


T(0)  „  _SL.<0,0)  u 

(36a) 

T<-">  .  „ 

(36b) 

2b 


This  stress  distribution  now  generates  a  forward  wave  of  amplitude 

-f-  .  -j. 

AA  and  a  backward  wave  of  amplitude  AA  (Fig.  7).  AA  is  In  quadrature  to 
A+,  leading  to  a  phase-shift  of  the  forward  wave;  the  phase-shift,  A0,  per 


periodic  section  is  given  by, 
A0  =  AA+/A+ 


(37a) 


and  the  reflection  coefficient,  r,  per  periodic  section  is  given  by, 
.  + 


r  =  AA  /A 

A0  and  r  are  related  to  the  stress  t  by  the  normal  mode  equations: 


.  itup  *  (0) 

a  0  =  u  •  t 


(37b) 

(12) 

(38a) 


r  =  4P2  u  *  C^"N^  (38b) 

a 

where  is  the  power  carrier  by  the  normal  mode.  Using  t^°^  and  t^ 
from  Equations  (36)  and  inserting  the  subscripts  running  over  x,y,z  ex¬ 
plicitly  (which  we  dropped  earlier  for  convenience),  we  have, 

z 

A0  =  -jp  ^  SL'<°’0)  Z  (39a) 

i,  j=x 


where 


r  =  -jp  ;  SL'f:N’0)  Z' 

Lj  lj 


\ 

L 

i,  j=x 


ij 


(39b) 


•k 

<iiu,  u. 

_ i _ i 


ij  4P 


(40a) 


a 


tuu.u. 
.  -  1  J 


ij  4P. 


(40b) 


Knowing  A0  and  r  we  can  deduce  AZ/Z  and  jB  for  the  transmission  line  model 
(Fig.  4).  By  cascading  the  transmission  matrices  of  successive  periodic 
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sections,  we  can  then  obtain  the  dispersion  relation  for  the  overall  array. 
This  normal  mode  approach  gives  the  same  results  as  the  direct  solution  of 
the  secular  equation  for  small  perturbations,  but  is  easier  to  implement. 
Moreover  it  yields  the  equivalent  transmission  line  parameters  directly. 

We  will  now  separate  the  first-order  and  the  stored  energy  effects 
in  Equations  (39)  using  the  expressions  for  SL'  given  in  Equations  (30). 


(1)  First-order  effect: 

The  first-order  effect  is  obtained  from  Equations  (39)  by  using 
only  the  first  term  (SL)  for  SL'  from  Equations  (30).  This  is  the  result 
obtained  from  the  coupled  mode  approximation.  In  this  approximation  we  have, 


A®  -  -j  (Pp)  (Ph)  )  lJ°’0)  z. . 

L  LJ  1J 


(41a) 


l,  j=x 


P 


r  =  - 


^  ( -N, 0) 

j  (Pp)  (ph)  )  L.  .  Z\  . 

L  iJ 

i»j=x  3 


(41b) 


Here  we  have  used  SL  in  place  of  SL'  in  Equations  (39)  and  used  the  fact 


that  (Equation  (22a)), 

£  SI.!0'01  ,  (Ph) 

4-  tj  j 

j j=x 


z 

\  ,(0,0) 

,  L.  .  u . 

/  ■  1J  J 


(42a) 


since  z 

L  Sij 
j-x 


(42b) 


u  being  the  particle  displacements  of  a  normal  mode  for  the  semi-infinite 

substrate.  In  Equation  (41b)  we  have  replaced  SL^'^’0'*  by  (Bh)  . 

il  ij 

since  S(;N’0)  =  0. 


i  J 


Equations  (41)  show  that  the  phase-shift  and  reflection  obtained 


from  the  coupled  mode  approximation  is  linear  in  (h/>). 


Using  the  expres- 


sions  for  L  from  Appendix  1  this  is  seen  to  be  the  same  result  as  in 
Reference  9. 


(2)  Stored  energy  effect: 

The  stored  energy  effect  is  obtained  from  Equations  (39)  by  using 
the  second  term  for  SL'  from  Equations  (30): 


2  V  H 


(0,0) 


U  =  i  Op)  (Sh)  )  ii  z. . 

j  — n —  1 1 


i,  j=x 


e 


where 


2  V  h<-N>°> 

=  J  <f»P>  (Bh)  l  Vl _  Z'ij 

i,j=x  P 


(0,0)  V  Y  T(°.n)  r(n,m) 

Hij  V  Lik  Gk  l  LAj 

n,m  X#k=x 

h(-N,0)  \  j ( -N, n)  (n,m)  L(m, 

lj  i.fex  ik  k* 


(43a) 


(43b) 


0) 


(44a) 


(44b) 


Hera  [SL]  has  been  replaced  by  (0h).[Lj  since  only  the  matrix  elements  be¬ 
tween  different  space  harmonics  are  involved,  for  which  [S]  =  0.  In  Equations 
(44)  the  summation  over  n,m  is  carried  out  over  as  many  space  harmonics  as 
desired. 

The  [L]  matrix  describing  the  layer  is  calculated  assuming  that 
it  is  a  thin  strip  uniformly  excited  across  its  thickness.  This  is  an  invalid 

l_l_ 

assumption  at  high  spatial  harmonics,  since  the  penetration  depth  at  the  n 
harmonic  is  of  order  ~  1/fV^.  Thus,  the  matrix  element  is  overestimated 

by  a  factor  jV^h  at  high  values  of  n.  An  exact  analysis  of  a  finite  thickness 


strip  is  rather  complicated;  a  good  approximation  is  to  view  the  strip  as  a 
single-mode  lossy  transmission  line  in  the  direction  perpendicular  to  the 
surface.  The  transmission  line  is  "shorted"  (zero  stress)  at  the  top  and  [ L ] 


represents  the  "impedance"  (stress/displacement)  at  a  distance  h  from  the 
load.  From  transmission  line  theory  the  impedance  Z  depends  on  h  through 
the  relation, 

z  =  z  exp(3(n)h)  -  exp(-p(n)h) 

°  expO^^h)  +  exp(-@^n^h) 

Z  3(n)h  for  g(n)h  «  1 
o 

where  is  some  constant  and  the  decay  constant  along  y  is  assumed  equal 
to  3^  for  the  nth  harmonic.  The  thin  strip  approximation  is  valid  when 
3*'n^h  «  1  and  yields  an  impedance  linear  in  h.  An  approximate  correction 
for  L^m,n^  £S  obtained  by  multiplying  it  with  the  factor, 

1  exp(3^n’>h)  -  exp(-g^h) 

P  h  exp( p  h)  +  exp( -3  h) 

This  correction  factor  is  clearly  not  precise;  however,  the  end 
result  is  quite  insensitive  to  the  precise  nature  of  the  correction.  This 
is  because  the  chief  contribution  to  the  stored  energy  effect  comes  from  the 
lower  space  harmonics  which  are  unaffected  by  the  factor;  the  higher  space 
harmonics  contribute  very  little  individually.  The  inclusion  of  the  factor 
does  serve  to  eliminate  spurious  contributions  arising  from  the  cumulative 
effect  of  a  large  number  of  space  harmonics. 
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G .  Results : 


(1)  First-order  effect  (Equations  (41)): 

Equations  (41)  give  the  first-order  A0  and  r  as  obtained  from 

the  coupled-mode  approximation.  Using  the  expressions  for  1.^’°^  and 

ij  ij 

from  Appendix  I  and  noting  that  ft  =  N.-^,  A0  and  r  can  be  evaluated  in  a 
straightforward  manner.  The  reflection  coefficient  is  equal  to  AZ/Z  for 
the  transmission  line  model.  The  results  agree  with  Reference  9  and  the 
calculations  are  not  repeated  here.  For  grooved  arrays,  A0  is  found  to  be 
zero  thus  predicting  no  shift  in  the  stopband  center.  For  electrodes  of  a 
material  different  from  the  substrate  a  non-zero  A0  is  obtained,  but  this 
is  of  order  ~  (h/X)  and  corresponds  to  a  change  in  velocity  of  the  wave 
under  the  electrode  material  equal  to  that  due  to  a  uniform  layer  of  the 
electrode  material.  Thus,  this  approximation  yields  a  repetitively  mis¬ 
matched  transmission  line  model  with  different  velocities  in  the  gaps  and 
under  the  electrodes.  The  change  in  wavenumber,  A ft,  under  the  electrodes 
is  obtained  from  A0  using, 


Aft  _  A0 


(45) 


where  T]  is  the  ratio  of  the  electrode  width,  a,  to  the  period,  p  (Fig.  5) 


(2)  Stored  energy  effect  (Equations  (43)): 

Equations  (43)  give  the  A0  and  r  due  to  the  stored  energy  effect. 

(i)  Phase  shift  at  the  fundamental: 

The  phase-shift  at  the  fundamental  was  calculated  for  aluminum 
electrodes  on  quartz.  We  compared  our  results  to  the  experimental  data  provided 
to  us  by  Rome  Air  Development  Center  for  a  double  electrode  transducer  on  ST-X 
quartz.  The  electrodes  were  made  of  80  K  of  chromium  and  820  X  of  aluminum. 
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The  line  width  was  1.1729pm  corresponding  to  a  fundamental  wavelength  of 

9.3832pm.  The  fundamental  frequency  of  a  double  electrode  transducer  is  at 

half  the  frequency  of  the  first  mechanical  stopband.  We  calculated  the 

change  in  velocity  of  the  surface  wave  in  propagating  under  the  electrodes. 

The  results  are  shown  in  Table  I.  The  piezoelectric  effects  on  ST-X  quartz 
2 

are  smaller  (K  / 2  =  .00055)  and  have  been  neglected.  It  is  seen  that  a 

Table  I:  A v/v  for  a  double  electrode  transducer  on  ST-X  quartz. 

Calculated 
,st 

1  order  Stored  Energy  Total 

.00065  .00385  .0045 

major  part  of  the  lowering  of  velocity  is  due  to  the  stored  energy  effect 
rather  than  the  linear  first-order  effect. 

figure  8  shows  the  results  for  grooves  in  ST-X  quartz  with  Ky 

given  by, 

H  =  (Pp)  .  Kv  .  (h/A)2  (46) 

The  experimental  points  are  taken  from  Reference  13  and  show  fair  agreement 

It  was  pointed  out  to  the  author  by  Melngailis  that  the  large  scatter  in 

the  experimental  data  could  possibly  be  due  to  varying  groove  shapes;  the 

present  analysis  assumes  a  rectangular  shape  thov'gh  other  shapes  can  be 

accounted  for  by  a  modification  in  1 1, ]  . 

The  phase-shift  Cor  grooves  in  Y -X  lithium  niobate  was  calculated 
2 

to  be  25(h/A)  at  h/A  -  .01.  This  is  in  agreement  with  the  measured  result 


Measured 


.0048 


in  Reference  14. 


(it)  Reflection  at  the  second  harmonic: 

The  reflection  per  strip  at  the  second  harmonic  for  grooves  in 

o 

Y-Z  lithium  niobate  was  found  to  be  22(hA)  at  hA  -  .01  which  is  about 
half  of  the  measured  value  in  Reference  6.  The  reason  for  this  discrepancy 
is  not  clear  to  us. 

The  reflection  coefficient,  r  per  strip  was  calculate!  at  the 

second -harmonic  stopband  for  aluminum  electrodes  on  ST-X  quartz  and  for 

grooves  in  ST-X  quartz.  The  results  are  shown  in  Fig.  9  where  K  is 

,  R 


Fig.  9.  vs  h/X  at  the  second  harmonic  for, 

(a)  Grooves  in  ST-cut  quartz 

(b)  Aluminum  electrodes  on  ST-cut  quartz. 
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given  by. 


r  =  KR  .  (h/X)2  (*7) 

This  result  has  not  been  compared  to  experimental  data. 

H.  Summary: 

A  theoretical  model  for  surface  wave  propagation  in  a  periodic 
array  of  thin  strip  overlays  has  been  developed.  The  "stored  energy"  effect 
is  modeled  by  extending  the  analysis  beyond  the  coupled  mode  approximation. 
The  results  are  discussed  for  aluminum  electrodes  and  grooves  in  ST-X  quartz 
and  Y-Z  lithium  niobate. 
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4.  PIEZOELECTRIC  SCATTERING: 

A.  Introduction: 

The  propagation  of  surface  waves  in  an  infinite  periodic  array 

of  thin  conducting  strips  has  previously  been  analysed*"^  ^  ;  and  an  equi- 

(19) 

valent  "mixed  circuit"  model  has  been  obtained  by  matching  dispersion 
relations  near  the  stop-band.  In  this  section  we  adopt  a  different  approach 
that  leads  directly  to  the  scatter  matrix  of  a  single  strip  in  a  periodic 
environment.  We  assume  (Fig.  10)  that  an  unperturbed  surface  wave  is  inci¬ 
dent  on  a  strip  within  a  periodic  array.  The  incident  wave  has  a  surface 
potential  0^  given  by, 

0g  =  0+  e'^Z  (48) 

where  p  is  the  wavenumber  of  the  unperturbed  wave.  The  charge  induced  in 
the  strip  by  this  surface  potential  is  determined  from  the  electrostatic 
equations.  This  charge  is  then  used  as  the  source  term  in  the  normal  mode 
equations to  determine  the  change  in  the  incident  wave,  A0+  and  the  gener¬ 
ated  reverse  wave,  A0  .  The  scatter  matrix  describing  a  strip  in  a  periodic 
array  is  thus  obtained.  Finite  arrays  (neglecting  and  effects)  are  then 
readily  analyzed  by  cascading  the  transmission  matrices  of  successive  strips. 
Since  the  formulation  is  on  a  strip  by  strip  basis,  assuming  local  periodicity, 
slow  variations  in  periodicity  and  metallization  ratio  may  be  accounted  for. 

The  advantages  of  this  approach  are: 

(1)  The  scatter  matrix  of  each  strip  is  determined  directly 
at  all  frequencies  so  that  the  model  is  valid  outside 
the  stopbands  as  well. 

(2)  Analytical  expressions  are  obtained  for  each  scatter 
matrix  element  as  a  function  of  metallization  ratio 


and  frequency. 
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(3)  The  analysis  can  also  be  extended  to  include  bulk  wave 
generation  by  using  the  charge  induced  on  the  strips  as 
the  source  term  in  the  normal  mode  equations  for  the  bulk 
waves . 

The  formulation  presented  ..ere  is  directly  applicable  to  periodic 
arrays  that  are  electrically  periodic  as  well.  This  includes  open-  and 
short-circuited  reflector  arrays  and  multistrip  couplers.  An  extension  of 
the  analysis  to  transducers  with  the  electrodes  arbitrarily  interconnected 
is  being  developed. 


Z  -  0  or  -  tin 


Fig.  10.  Reflection  and  transmission  of  surface  acoustic  waves  at  an 
electrode  in  a  periodic  array. 


B.  Charge  Distribution  Induced  In  a  Periodic  Array  by  a  Propagating 
Surface  Waver 

In  this  section  we  will  determine  the  Fourier  components  of  the 
induced  electrostatic  field  in  a  periodic  array  due  to  an  incident  surface 
wave.  For  convenience  dimensionless  variables  s  ,  s  and  9  are  introduced. 
2ttz 


9 

P 

s  .  ^ 
S  2tT 

s  '«  £E 
S  2tt 


(49a) 

(49b) 

(49c) 


where,  3  is  the  wavenumber  of  the  surface  wave 
and  9  =  k  +  N.  2tt/p( 

N  being  a  positive  integer  such  that  0  <  k  <  2xr/p.  It  follows 

that , 

s  =  s  +  N  (49d) 

where  0  <  s  <  1 . 

The  surface  potential  due  to  the  unperturbed  wave  is  given  by 
-j  (s-t-N )  9 


0  =  <4  e 

s  + 


(50a) 


In  the  presence  of  the  periodic  array  of  conducting  strips  an  additional 
electrostatic  field  is  generated.  This  field  is  periodic  in  z  and  its 


surface  potential  <b  is  written  in  a  Fourier  series  as, 


V9)  =  L  K  e 

rt*  -oo 


-j  (s  +  N)  9 


(50b) 


The  total  surface  potential  is  the  sum  of  the  wave  potential  and  the 
induced  electrostatic  potential. 


«(8)  -  »  9  +  }  »  9 
+  nitm  n 


(51) 


38 


The  tangential  electric  field,  at  the  surface  is  obtained  from  (C>1) 


E ,<«)  -  Jin  .--Is  9 

Z  p 


(st-N)®  e  +  '  (s+n)0  e  J 

+ 

n=  -<*> 


-  in9 


(  52a  > 


The  electrostatic  part  of  the  field  a  Iso  gives  rise  to  a  norma  1  electrical  dlsp  lace- 
men  t,  Dy  at  the  surface  (or  a  charge  distribution)  that  is  written  from  Laplace  1  s  equation 


D  (9)  =—(€+€  )e  >  s  (stn)$  e  J“  WJ 

”  P  P  o  L  _Lj  _  n  n  J 


-jn  0"; 


( 52b) 


n=  -<*> 


where,  S  *  sign  of  (s+N) 
n 

=  +1  for  n  >  0 

=  -1  for  n  <  0 

since  0  <  s  <  1. 

and  is  the  effective  permittivity  of  the  substrate  and  £q  is  the  permit¬ 
tivity  of  the  medium  above  the  substrate. 

The  Fourier  coefficients  0  have  to  be  determined  such  that  the 
boundary  conditions  imposed  by  the  array  are  satisfied  (Figuren): 

let  <  a 


Ez(9)  =  0 
Dy(9)  =  0 


A  <  0  <  TI 


(53a) 

(53b) 


where  A  =  Hit 

T]  being  the  metallisation  ratio. 

These  conditions  are  written  from  Equation  (52a)  and  (52b)  as. 
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where , 
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(s+N)T 


(54c) 


39 


(/  j<S,NH> 


A 0  o>(s  ,  N,°- 


0-  A 


0-  i  A 


-  » 


(</>,  +  Af/>(  )e 
A  =  7f 


I  n 


0--  0 


|(S  IN)0 


11.  Same  as  Figure  10  with  dimensionless  variables. 


It  is  shown  in  Appendix  IV  using  the  properties  of  the  Legendre  polynomial 
that  conditions  (54a)and  (54b)  are  satisfied  if  we  write, 


N 

A  =  )  or  S  P  .(cosA)  +  or,  ,  S  P  (cosA) 

n  m  n  rtf-m-N-1  N+l  n  n 


(55) 


where  or ,  ,  is  any  arbitrary  number  and  to  or.  are  determined  from 
N+l  1  J  ON 

the  following  equations: 


a0  =  "  2 


(56a) 


X-l 


at,  =  -  /  at  P.  (cosA)  ,  £  =  1  to  N 
i  t—  m  i- m 
m=0 


th 


(56b) 


Here  P^  represents  the  n  Legendre  polynomial  and  for  negative  n  is 

obtained  from  the  relation  P  ,  =  P  .  The  charge  distribution  in  the 

-n-1  n 

array  is  written  from  Equations  (52)  and  (55), 

Nf  1 
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)  /Of  )  P  ..  .  (cosA)e  J 
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m=0  n*  •“ 

(57) 


0 


The  terminal  voltage  and  current  at  the  strip  are  shown  in  Appendix  V  to  be, 

jif  i 

N-nH-1  ti  p  (-cosA) 

V  =  (sfN)0,  )  (-1)  Cf  - -  (58a) 
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J  v  p  o'  +  m  N-m+sv 


(58b) 


where  W  is  the  length  of  a  strip  perpendicular  to  the  direction  of  wave 

propagation.  aQ  to  af^  are  uniquely  known  from  Equations  (56a)  and  (56b) 

has  to  be  chosen  so  that  the  correct  ratio  V/l  =  Z  is  obtained  from 

Equations  (58a)  and  (58b)  . 

It  is  convenient  to  divide  of  into  two  parts: 

£*fl 


=  aN+l  + 


a 


N+l 


(59) 
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where  a  is  chosen  such  thac  C.ne  scrip  v outage  ts  zero.  this  correspor.cs 
co  che  case  when  all  che  scrips  are  shorced  together.  ,-rom  equation  < 5ba  >  , 
-  N-m  P„  ( -cos  A) 

a '  =  )  (-1)  m  N'CH'S  (60) 

^  _ 

?  ( -cos  A) 

S  -  1 
V 

i  is  also  obcained  crotn  Equation  (68a): 


.  . V  Vs  inns 

(s-r.  )  aM^l~  n  ?  (-COS  A I 

3  -1 


Pot  I  S  (I  lf*(  lri< dl ) 
V 


0-  -  IT  0  -  t  V 


Fig.  12.  Scattering  at  an  electrode  in  a  periodic  array. 
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U.  Evaluation  of  Scatter  Matrix  for  a  Single  Strip: 

In  this  section  we  will  use  the  induced  electrostatic  field  de¬ 
rived  in  equations  (56) through  (61)  to  obcain  the  scatter  matrix  of  a  single 
Strip.  To  start  with  we  define  a  (3x3)  matrix  A  that  relates  the  fields  at 
the  two  acoustic  and  one  electrical  ports  (Fig.  12).  All  values  are  refer¬ 
enced  to  the  center  of  che  electrode. 


(62a) 


The  (2x2)  scatter  matrix,  S  relating  the  two  acoustic  ports  is  easily  obtained 
from  the  A-matrix  by  setting 


I  =>  V/Z  (62b) 

where  2  is  the  external  load  at  each  strip.  Using  equation  (62b)  in  equation 
(62a)  we  get, 


(63) 


where , 


S  a  A  til  tlL. 

u  '  ii  a,.,-  n'z') 


3  3 


-  ,  A13  A32 

12  i2  ‘  A„-  (1/2  ) 


S.,  =  A,.  - 


33 

A23  A31 


a33-  (I1:' 


3,,  3  A .  - 


A23  a3; 


'33 


( 64a) 


(64b) 


( 64c) 


(64d> 
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We  will  now  evaluate  the  elements  of  the  A-matrix.  The  induced 


electrostatic  field  was  derived  in  the  last  section  assuming  an  incident  wave 
only  from  the  left;  that  is,  $  was  set  equal  to  zero.  We  will  thus  evaluate 
the  A-matrix  elements  in  the  first  and  third  columns  (Equation  (62a)).  The 
elements  in  the  second  column  are  written  from  the  symmetry  of  the  acoustic 
ports: 

■X  =  A 
12  21 

A0O  =  Au  (65) 

A32  =  A31 


Using  equations  (58b),  (59),  (60)  and  (61)  the  terminal  current,  I 
is  written  as  the  sum  of  two  parts--one  proportional  to  0  and  the  other  pro¬ 
portional  to  V.  Comparing  with  equation  (62a)  ,  we  obtain  and  A^. 


_  N 

A31=  jaJ(V  €0)W-2TJ(S+N)-lnam'vPN-n^s(COsA) 

m=0 


(-‘/-"^(cosA)  ^s(-cosA>- 

P  . (-cosA)  /_ 

s-1 


A 


33 


2J«(V 


€q)W  .  sinTTs 


P  . (cosA) 

s-1 _ 

P  . ( -cos  A) 
s-1 


(66) 

(67) 


A  ^5  represents  the  strip  current  with  V  =  0,  that  is  with  the  array  shorted 
while  rePresents  the  strip  current  due  to  a  voltage  V  with  no  incident 

wave  3  .  It  will  be  noted  that  since  the  fields  were  assumed  periodic,  a 
voltage  V  means  that  the  voltage  at  each  strip  has  the  same  magnitude  V  but 
is  delayed  in  phase  with  respect  to  the  previous  one  by  2stt. 

Our  next  step  is  to  evaluate  A^,  A^>  A^j  and  A?^-  This  is  done 
by  using  the  induced  charge  distribution  0^,(9)  obtained  in  the  previous  section 
(Equation  (57) to  determine  the  change  in  the  forward  wave,  A®^  and  the  gener¬ 
ated  reverse  wave,  A3  from  the  normal  mode  equations:^1* 
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dK 

=  iE 

1  j  Av  j 

d3 

2n 

€P+  So  I  v  1 

d® 

-  -iE 

1  J  Av 

IT 

2n 

6  +*  €n  v 
p  0  ‘ 

-j  (s+N)  3 


(68a) 


(08b) 


Integrating  over  one  period  trora  -n  to  +n  and  assuming  constant  over 
the  period, 

a  tf+1 

!.'v!  “  *  -Arn-N)9 

— —  =  jl — )  (s+N)-;  a  ;d9  ?  .(cosA)  e  J 

a  v  L~  m  j  n^m-M-L 

+  1  tn=D  -n  n=  -® 
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=  jTT 


w-i 
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since 
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d9  e  ji'n  »  2ti  if  n=N 


“  -rr 

A<fi_ 

T" 


=  0  otherwise 
N+l 

*TT 
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L, 
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=  j  j~ I  ,2tt(s+N) 
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\ 
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v  a  ?  (cosA) 

ai  2.V+2s  -m 


(69b) 


since  !  d9  )  ?  ,  (casA)e  a  2tt  .  Pnu_9c  (cosA) 

u  _  l-j  cmn-N-l  zN+zs-m 

-TT  n=  -  = 


as  shown  in  Appendix  VI.  Using  Equations  (59),  (60)  and-  (61)  to 
a  ,  in  equations  (69a)  and  (uvb).  A®  and  A®  are  written  as 

iN  i  •A»  - 

two  parts--one  proportional  to  ®  and  the  ether  proportional  to 
yields  the  desired  A-matrix  elements. 


replace 
the  sum  o 
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A  =  1  +  j  I—  I  .  2tt(34-N) 

1  1  !  V  1 


ct  '  P  .  (cos A)  + 
— ~  m  -  m-l 
m=u 
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( 7  Od ) 


The  matrix  elements  derived  earlier  are  reproduced  here  for  convenience 


A  =  itu  (€  €j  W.  2tt(s-*-M)  . 
31  ?  0 
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( 70e) 


?  .  ( cosA) 

A  =  2j«i(€  +■  €  )  V  sinrrs  -^Li- - - 

?  0  ?  1  (-cos  A) 


(  70  f ) 
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A21 

(70g) 
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(70h) 

A3L 
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Equations  (70)  give  the  complete  scatter  matrix  for  a  strip  in  a  periodic- 


array  once  the  a'S  are  determined  from  Equation  (56).  A  few  of  the  a'S  are 
m  m 


calculated  here  for  convenience. 
1 

a0  ~  '  2 


(71a) 


=  J  cosA 


1  .  2.  .. 
a2  =  Z  (C0S  A_1) 


1  3 

=  -  (cos  A  -cosA) 


(71b) 


(71c) 


(  7  Id ) 


Equations  (64)  are  used  to  obtain  the  scatter  matrix  of  a  single  strip  for 
a  given  load  termination  Z,  from  the  A-matrix  derived  above. 


D.  Examples: 

We  now  have  analytical  expressions  for  the  scatter  matrix  elements 
for  all  frequencies  and  metallization  ratios.  In  Ref.  19  the  reflection 
coefficient  (S^)  pec  strip  was  computed  at  various  stopbands  for  open-  and 
short-circuited  arrays.  In  this  section  we  will  compare  our  results  with 
Ref.  19  at  the  fundamental  and  second  harmonic  stopbands. 

(a)  Fundamental  stop-band:  (s  =  .5,  N  =  0) 

At  the  fundamental  stop-band  s  =  .5,  N  =  0 
Equation  (71a)  in  Equations  (70)  we  get  the  relevant  A-matrix  elements: 


Using  c*Q  =  -  1/2  from 


11 


.  Av 


T.  r  P<5(-cosA) 

•  2  EosA  +  p"  (-cos'A)j 


-.5' 


A13  =  2j 


Av _ 

v  [  P  .(-cosA) 


A31  -j«(6  +  60)W.f  [p>5(cosA)  + 


P  (cosA) 

A„  *  2jm(€  +  €n)W.  •= — '■ — - 7T- 

33  J  p  0'  P  ^(-cosA) 


P  (cosA)P  (-cosA)- 

—  •  -J  •  D 

P  ^(-cosA) 


(72a) 
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(72d) 
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Here,  che  Legendre  polynomials  P^cosA)  =■  cosA  and  PQ(cosA)  =  1  have  been 
used.  From  Equations  (64)  we  now  write  the  reflection  coefficient  for  the 
open-  and  short-circui ted  array  by  setting  Z=“  and  Z=0  respectively. 
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(73b) 


Equations  (73)  provide  the  reflection  coefficient  of  an  open-  and  short- 
circuited  array  at  the  fundamental  stop-band  as  a  function  of  metallization 
ratio.  These  agree  with  che  plots  in  Ref.  19 

(b)  Second  Harmonic  stop-band:  (s=l,  N=0) 

For  che  second  harmonic  stop-band  we  may  use  s=0,  N=  1  or 
s=l,  N=0;  either  choice  yields  the  same  result.  In  this  case  A^=  0  so  that 
5^1  is  independent  of  the  load  Z  which  is  expected  intuitively.  We  have  for 
che  second  harmonic  stop-band, 


S11  =  A11 


.  Av  TT  .2. 
=  J  — .m  sin  A 
v  2 


(74) 


for  all  Z . 
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E.  Summary: 

The  scatter  macrix  of  a  single  scrip  in  a  periodic  array  has  been 
derived  analytically  for  all  frequencies  and  metallization  ratios.  Analytical 
expressions  for  the  reflection  coefficient  of  a  strip  in  an  open-  and  in  a 
short-circuited  array  at  the  fundamental  and  second  harmonic  stop-bands  are 
derived  as  special  cases  of  the  general  results. 
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5.  ELEMENT  FACTOR  FOR  PERIODIC  TRANSDUCER: 

A.  Introduction: 

A  large  number  of  practical  transducers  are  modeled  as  a  periodic 

array  of  metallic  strips  with  arbitrary  impressed  voltages.  It  has  been 

(21  22) 

shown  ’  that  the  response  of  such  periodic  transducer  is  written  as  the 
product  of  an  "element  factor"  and  an  "array  factor".  The  element  factor 
was  obtained  in  Ref.  21  using  the  numerical  tables  of  Ref.  23  which  take  the 
effects  of  two  neighboring  electrodes  into  account.  In  this  paper  analytical 
expressions  for  the  element  factor  are  obtained  for  all  frequencies  and 
metallization  ratios  with  the  effects  of  all  neighboring  electrodes  accounted 
for.  The  response  of  the  periodic  transducer  (neglecting  regeneration  and 
reflection)  is  thus  obtained  analytically  for  all  frequencies. 

We  first  consider  a  phased  array  transducer  (Fig.  13)  in  which  the 
voltages  on  all  the  electrodes  have  the  same  magnitude  V  but  the  phase 

Ve^2koP  VeikoP  V  Ve'j^P  Ve~i2koP 

- 


KP-1705 


ig.  13.  A  phased  array  transducer  with  a  uniform  phase  progression  of 

k  p  per  period. 

o 

progresses  uniformly  along  the  array  at  the  rate  of  kQp  per  period  where 

p  is  the  period.  The  voltage  on  the  n^  electrode,  V  is  written  as, 

n 

V  -  V  exp( -  ink  p) 
n  or 


-  V  exp(-in2tTs  ) 
o 


(75) 
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(76) 


where  s  is  a  dimensionless  variable  related  to  k  by, 
o  o 

s  =  k  p/2n 
o  o 

It  will  be  noticed  from  Equation  (75)  that  adding  any  integer  to  s  yields 

the  same  electrode  voltages  so  that  we  may  without  loss  of  generality 

restrict  the  values  of  s  such  that 

o 

0  <  s  <  1  (77) 

o 

(24) 

The  charge  distribution  in  this  transducer  is  known  to  be 


Q( 9)  =  —  (6  +  €  ) 
v  P  P  o 


s  inns 

o 


(-cos'A) 

o 


P  (cosA) 
n 


-j (s  +  n) e 
e  o 


(78) 


where  (1)  0  is  a  dimensionless  variable  related  to  z  by 

0  =  2nz 
P 

(2)  6  is  the  effective  permittivity  of  the  substrate  and  6 
v  p  o 

is  the  permittivity  of  the  medium  above  the  surface. 


(3)  A  =  7]n 

T)  being  the  metallization  ratio  and 

(4)  is  the  Legendre  function  of  order  v  (integer  or  non- integer ) . 
From  Equation  ( 78)  we  note  that  the  Fourier  transform  Q(s)  of  Q(0)  has  delta 
function  components  at  all  s  given  by  s  =  sq+  n,  n  being  an  integer. 

Q(s)  =  Q(sq)  .  Pn(cosA)  (79a) 


where 

and 

Here  s  is 


s  =  s  +  n 
o 


S'V  '  f  <«p  +  V  r 


SI  HITS 


S  -  1 
o 


( -cos A) 


(79b) 


the  usual  spatial  harmonic  k  scaled  by  p/2ir  to  make  it  dimen¬ 


sionless.  Fig.  34  shows  a  plot  of  Q(s)  against  s  for  the  phased  array 
transducer  with  electrode  voltages  given  by  (7r)). 
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Fig.  1m-.  1  ourier  transform  of  charge  distribution  for  a  phased  array 

transducer. 

The  electrode  current  at  each  strip  has  the  same  magnitude  but 

with  a  phase  progression  along  the  array.  The  current,  D  at  the  nth 

n 

strip  is  written  as 

Dn  =  Do  exp('js02TTn)  (80a) 

where:  is  obtained  by  integrating  ju>  Q(0)  over  the  width  of  a  strip 
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jo;  W 
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+A 

i  d0  Q(0) 
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(6  +  €  )  .  2V 
P  o' 
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ps  _ 2 (cosA) 

o _ 

o  Pg  _1(-cosA) 


(80b) 


o 

W  is  the  length  of  a  strip  in  the  direction  perpendicular  to  the 
wave  propagation. 

Mow,  we  wish  to  determine  the  charge  distribution  in  a  single 
tap  transducer  (Fig.  15)  since  this  yields  the  element  fac tor . ^2 1 ’ 22^  jt  £s 
easily  seen  that  this  is  obtained  by  superposing  the  charge  distributions 
in  an  infinite  number  of  phased  array  transducers  each  with  V  =  1,  but  with 
,;(>  '■’arving  uniformly  I  rmn  0  to  1.  The  voltage  on  the  n1*1  electrode  is  then 
obtained  from  liquation  (75)  as 
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IV 


Fig.  15.  Electrode  currents  in  a  single  tap  transducer. 

1 

V  i  ds  exp(-jn2ns  ) 

n  J  o  J  o 7 

o 

=  1  for  n  =  0 
=  0  for  n  ^  0 

which  clearly  describes  the  single  tap  structure  in  Fig.  15. 

Hence  the  Fourier  transform  of  the  charge  distribution  in  a 
single  tap  structure  is  a  continuous  function  of  s  given  by, 

2(6  +  6  )  sintts 

ct(s,A)  =  - ^  _i(.c°sA)  •  Pn(cosA)  (81) 

o 


where  s  =  s  +  n 
o 

1/2  — 

such  that  0  <  sq  <  1  and  n  is  an  integer.  Fig.  16  shows  s  a  (s,A) 
plotted  against  2s,  for  different  values  of  the  metallization  ratio.  These 
agree  reasonably  with  the  numerical  plots  in  Fig.  15  of  Ref.  21  where  2s 
corresponds  to  f/f^.  The  discrepancies  are  possible  because  only  two  near¬ 
est  neighbors  were  accounted  for  in  the  numerical  method. 

The  capacitance  of  periodic  transducers  was  obtained  in  Ref.  21  in 

terms  of  the  electrode  currents,  D  (Fig.  15)  in  the  single  tap  transducer. 

n 

These  are  readily  obtained  from  F.quation  (80)  by  integrating  over  8^  from 
0  to  1. 
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In  general  the  integral  has  to  be  obtained  numerically. 


Bui  for 
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cally  as, 
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then  obtained  analyti¬ 


cal 
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(83b) 


It  is  seen  that  while  is  positive,  all  other  D^'s  are  negative  as  expected. 
It  may  be  shown  that,  as  expected,  the  sum  of  all  the  electrode  currents  is 
zero .  os 


\ 


n=-« 


D 

n 
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B.  Summary: 

In  Ref.  21  the  conductance  and  capacitance  of  a  periodic  transducer 
with  arbitrary  voltages  were  related  to  the  charge  distribution  and  the  elec¬ 
trode  currents  in  a  single  tap  transducer.  However,  these  quantities  were 
obtained  from  numerical  tables  that  take  only  two  nearest  neighbors  into 
account.  In  this  section  the  charge  distribution  and  the  electrode  currents 
in  a  single  tap  transducer  are  obtained  analytically  considering  the  effects 
of  all  the  nearest  neighbors,  making  the  analysis  more  accurate  and  more 
convenient.  The  overall  charge  distribution  in  a  periodic  transducer  is 
given  by  Equation  (5)  of  Ref.  21  as  (lie  product  (in  k-space)  or  the  convo¬ 
lution  (in  x-space)  of  the  voltage  distribution  and  basic  charge  distribution 
function  (BCDF);  the  BCDF  is  now  known  exactly  from  Equation  (81). 
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APPENDIX  I  :  DELTERV  and  DEI, TER/  1  rom  Scattering  I’.iram.- te rs  : 

In  this  report  the  scattering  parameters  have  been  derived  lor  a 
single  section  in  an  infinite  array.  In  this  appendix,  the  impedance  mis¬ 
match  AX  and  the  fractional  change  in  velocity  Av/v  for  an  equivalent  trans¬ 
mission  model  are  derived  for  the  S-parameters . 

Let  u>  be  the  radian  Irequeucy,  v  the  unperturbed  surface  wave- 

velocity,  l  the  length  of  one  section  and  Tj  the  metallization  ratio.  In  the 

equivalent  transmission  line  model,  the  metallized  sections  have  an  impedance 
and  velocity  different  from  the  unmetallized  sections  (Fig.  Al).  In  addition, 
two  susceptances  jB  are  included  at  the  edges  ot  the  metallized  section  in 
order  to  account  for  the  even  harmonic  reflection  which  is  not  predicted  by 
the  simple  mismatched  transmission  line  model. 

A  straightforward  analysis  shows  that, 

S  ej9  =  j  I AX  sinTjO  -  B  cosT]0j  (Ala) 

S12ej  9  =  l  -  jT,9  —  (Alb) 

where  9  =  uu l/v  and  AZ ,  B  are  normalized  to  the  characteristic  imped- 

o 

ance.  Using  these  equations  AX  and  Av/v  are  calculated  1  rom  and  It 

will  be  noted  that  at  the  fundamental  with  T)  -  .5,  cosTj0  =  0;  so  that  B  does 
not  enter  the  picture.  On  the  other  hand,  at  the  second  harmonic  sin"0  -  () 
so  that  AX  does  not  enter  the  picture.  At  intermediate  frequencies  or  for 
”1  other  than  .5  there  is  no  unique  solution  lor  AX  and  B;  however,  it  is 
possible  that  a  particular  choice  mighL  result  in  values  that  are  nearly 
constant  with  Tj  or  with  frequency. 

Av  AX.  ( 1 ) 

DEL'IERV  and  DKLTERZ  are  obtained  ;  rom  —  and  —  from  the  relations: 

v  X 

•> 

Lr 

DELTERV  -  A-  +  —  (Ale) 

2  v 

BELTER/  -  -  V  -  x  (Aid) 


Al 


9  =  oj-P/\j0 

rj  -■  Metallization  Ratio 
(v  ~  Radian  Frequeny 
v0=  Unperturbed  Surface 
Wave  Velocity 

KP  -  1853 

K'luiva  lc:iit  transmission  line  model. 

f  (z) 


z  -  -  p/2  z  r  d  p/2 


i’unci  ion  describing  geometry  of  the  strips. 


APPENDIX  LI: 


In  this  appendix  we  will  evaluate  the  stress-displacement  relation 
for  a  thin  periodic  array  of  strips.  I.et  us  first  consider  a  particle,  dis¬ 
placement  vector  with  a  single  spatial  frequency  g^n'. 

Lh  =  Il[n)  expl-jB(n)Zj  (A 2 a ) 

where  (3^n^  =  g  +  n  .  2n/p 


If  the  thin  layer  were  uniform  rather  than  periodic  the  stresses  t  .  would  also 


have  the  same  spatial  frequency  g 


(n) 


The  relation  between  t  and  U  can  then 


(n) 


be  expressed  by  a  constant  matrix  [A  ']  such  that 


rfn)  =  !  A(n)  U(n) 

i  .  1 1  i 

j  =x  J 


r  -  x,  y,  7. 


(A2b) 


The  matrix  elements  A^"^  are  known  from  a  first-order  perturbation  ana  lys is . ^ ^ 


A(n)  =  g(n) 

(S3(n)h) 
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(Vn)h) 

a' 

XZ 

xz 
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A(n)  =  g(n) 
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1  a* 
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A(n)  =  g(n) 
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(P(n)h) 

l-p’ 

other  A^^  are 

zero,  where 

,  •  ,.,2  /  l2  I 


(A3  a) 

(A3b) 

(A3c) 

(A3dl 

(A3el 


ij 

and  a'  ,  a'  ,  or'  ,  a'  have  the  dimensions  of  stiffness  and  are  given  by, 
XX  xz  zx  zz  J  ’ 


a’  =  S^/D 

XX  J 
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xz  5 


(A4M 


a; 
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(A4e) 

(A4f) 


[S'J  is  the  compliance  tensor  of  the  strip  material. 

Now  if  the  array  of  strips  is  periodic,  it  has  been  shown1, 
that  the.  stress  may  be  written  in  terms  of  the  particle  displacements  as, 


(3,4) 


where 
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(A6c) 


(A6d) 


(Abe) 


wn  in  Fig.  A2  .  The  terms  involving  ^ 
in  Initiations  (A6)  represent  delta  functions  at  the  edges  of  the  strip;  the 
physical  origin  of  these  delta  functions  has  been  explained  in  References 
3  and  4.  Multiplying.  Kquat  ion  (A4)  by  q  j  P  rn^  j,  ^  integrating  over  one  period 


and  using  text  liquations  (lb)  and  3a)  we  have. 


r 


+  p/2 


~  |  (z)  cxpl  j  (0""'  -  0'“')  zj  dz 


,(m)  „(n), 


(0h)  .  L(n,m)  =  - 

1J  P  -p/ 2  1J 

Using  (z)  from  Equations  (Af>)  ,  wu  get  L^n,m^ : 

L("’m)  =  0  l.p-c2  +  «•  (1  +  Sl^Hl  +  |  p 

xx  L  s  xx  0p  pp  J  m-n 

L(n,m)  =  p  cr*  (1  +  -2^)(1  +  (n-m) .—)  f 
xz  xz  Pp  Pp  m-n 


(A7) 


(A8a) 


(A8b) 
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(A8c) 
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(A8d) 
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(A8e) 


Where, 

+p/2 

Fm  n  =  n  1  exP^J 

m-n  p  Jp/2 


P(n))zj  dz 


Ti  stn  (m  -  n)  ttT| 
(m  -  n)  Till 


(A8f) 


T]  is  the  duty  factor  =  a/p,  p  is  the  surface  wave  wavenumber  and  c  is 

s 

the  surface  wave  velocity.  The  other  L^n,m')  are  zero. 

ij 


APPENDIX  III: 

In  this  appendix  we  will  evaluate  the  stress-displacement  relation 
for  a  semi-infinite  substrate.  As  discussed  in  the  paper,  the  substrate 
does  not  introduce  any  coupling  between  different  spatial  harmonics  (Equation 
(5)  of  the  text) . 

Let  us  consider  a  particle  displacement  vector  with  a  single  spatial 
frequency  =  0  +  n.2Ti/p.  If  we  assume  that  the  particle  displacements 

decay  into  the  bulk  from  their  surface  values  given  by  (A9)  as  an  exponential 


A5 


o  then  from  the  field  equations  we  get  three  values  of  a.  These 

correspond  to  three  elastic  modes  of  the  solid.  The  field  equation  is 


written  as 


Sc  3  Uk 

L  ijki  - - r- 


1 


Sx .  Sx 

J  l 


2 

p<u  U. 
1 


where,  C  is  the  stiffness  tensor  and  p,  the  mass  density  of  the  substrate. 
The  coordinates  x^,  x^,  x^  are  used  in  place  of  x,  y,  z  for  convenience.  If 
we  replace  the  derivative  operators  by, 

'VoXj  -  0 

o/dx^  -  a 
b/'fa.y  — 

liquation  (A9)  becomes  an  algebraic  equation.  The  determinant  is  a  sixth 
order  polynomial  in  a,  so  that  by  setting  it  equal  to  zero  we  get  six  values 

of  a  of  these  three  represent  exponentially  growing  solutions  and  are  re¬ 
jected.  Corresponding  to  each  of  the  other  three  we  get 

particle  displacement  vectors  whose  surface  values  are  (u^^  },  {u^  }  and 

[u(3)}.  The  total  particle  displacement  {u}  at  the  surface  is  a  linear 
combination  of  these. 


[ll]  |uj  .  (A} 


(A10) 


where,  I  u J  is  the  matrix  formed  by  arranging  [u^^},  and  {u^^}  in 
columns.  £a }  is  a  vector  whose  components  A^,  A^ ,  A^  are  constants  repre¬ 
senting  the  linear  combination. 

Now  corresponding  to  each  solution  [u^^},  [u^^},  [u^"^},  we  may 
calculate  the  stress  vector  at  the  surface  {T^^},  the 


total  stress,  t  is  written  as  a  linear  combination  ot  these. 


{t}  =  (Tj  [A]  (All) 

Combining  (A10)  and  (All), 

[t]  =  IT]  fu]"1  [U]  f A  1 2 ) 

comparing  with  text  Equation  (3b). 

(S ]  =  [T]  [u]'1  (A13) 


The  procedure  outlined  here  may  be  used  for  any  anisotropic  substrate  to 
obtain  the  matrix  [S^n’n^j  relating  the  stress  and  displacement  at  a  given 
spatial  frequency  .  For  an  isotropic  substrate  it  is  possible  to  obtain 

analytical  expressions  for  the  matrix  elements.  The  three  values  of  a  are 
obtained  as. 


or 


=  ^  =  ([P(n)]  *  k£) 


1/2 


(A  14a) 
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(A 14b) 
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(A14c) 


where, 
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L  \  +  2p 


(A15a) 


k? 

S  p 


(A15b) 


where  X  and  p  are  the  Lame'  constants  of  the  substrate.  The  matrix  |u]  is 
given  by. 
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The  [Tj  matrix  is  obtained  as. 
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From  (Alb)  and  ( A 1 7 )  we  get  the  [S]-matrix. 
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(A18d) 
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The  other  are  zero.  For  an  arbitrary  anisotropic  substrate  analytical 

expressions  are  not  obtainable,  but  the  matrix  elements  are  numerically 
evaluated  for  a  given  using  the  procedure  outlined  above. 


A8 


APPENDIX  IV: 


In  this  appendix  the  values  of  a  are  derived  so  as  to  satisfy 

m 


equations  (7)  and  (8) : 


e'jN9  +  )  a  )  S  P  „  ,(cosA)e'jn°  =0  lei  <  A  1A19) 
m  L>  n  n+m-N-1  *  1 

tn=U  n=  - 00 


IT  T  i. 

m=u  n=  - c 


Pn+m-N-l(COSA)e'Jn9  =0  A  <  |9I  <n 


(A20) 


It  is  known  from  the  properties  of  the  Legendre  polynomial  that1, 


)  S  P  (cosA)e  J  =0 

L,  n-ra  n-m 
n=  -cd 

w  . 

)  P  (cosA)e  ■'n  =  0 

L->  n-m 
n=  -ro 


0  <  A  (A21) 


A  <  0  <  u 


(A22) 


where  m  is  any  integer.  From  (A4)  it  is  seen  that  (A2)  is  satisfied  for 

any  choice  of  a’s. 

m 

In  order  to  satisfy  (Al)  we  choose  a’s  such  that, 

J  '  m 
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From  equation  (A3)  we  note  that  the  right  hand  side  is  0  for  [ Q ]  <  A. 


Equation  (A5)  yields, 
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where, 
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Since  3^  +  ^  -  0  it  is  apparent  from  Equation  (A6)  that  a^ 

arbitrary  value.  Leaving  this  term  out, 
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Equation  each  po..'r  (e  J  )  to  0  starting  from  e  we  get, 
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Equation  (A7i  is  written  as, 
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noting  that  P  .(cosAj  -  1  and  P  =  P. 
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APPKNIHX  V; 

In  this  appendix  the  terminal  voltage  and  current  at  a  strip  are 
determined  from  the  tangential  electric  field  and  the  normal  electrical  dis¬ 
placement.  From  Equations  (5a),  (7c)  and  (8), 
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A  10 
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where  Equation  (A5)  has  been  used. 


Integrating  (A9) , 
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Since  the  potential  is  uniform  over  a  scrip,  the  strip  voltage  V  is 


equal  to  0(0) . 
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using  the  derivation  in  Ref.  7.  The  strip  current  is  obtained  by  integrating 


Dy(8)  over  one  strip  width.  From  Equation  (10), 
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again  using  the  result  in  Ref.  7. 


Nfl 

,+A  “ 

/  a 

j  ) 

Ly.  m 

J  .  U 

m=U 

-A  n=  - 

Nt-l 

)  a 

.  2n  P 

m4t)  m 

C 

(cos  A) e 


■  j (s+n) 6 


(A30) 


APPENDIX  VI: 

In  this  appendix  the  following  integral  used  in  Equation  (20) 
is  computed: 

+TT  jb 


J  =  f  d9  \  P  (cosA)e'^  (2s+n+N)  9 

c.  n+m-N-1 


-tt  n=  -® 


,,  ,  T  ,  ..  2sm  n  (2s+n+N) 

We  have,  J  =  /  P  ,,  .(cosA).  — — - - 

’  m  n+m-N-1  (2s-*-n+N) 


All 


=  2  sin  2stt  .  ^  (-1) 

n=~  -°° 


n+2N-m+l  P  (cosA) 
n 


2s+n+2N-rM-l 


2  N-m-f  1  \  n  *  T  1  1 

=  2sin  2sn  .  (-1)  .  /(-I)  P^osA) ^Zs+ZN-nH-l  "  f^2s-n 

n=0 


mj 


2N-m  n 

2  jin  2stt  .  (-L)  P  .  (cosA)  .  — - - - 

2fH-2s-m  srnri(2N+2s-in) 


using  che  definition  of  Legendre  functions  P  for  non-integer  v. 
This  yields, 

J  =  2tt  .  P...  .  (cosA) 

2!$+-2s-m' 
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